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Abstract 


Composites offer unique opportunities in aerospace design. Apart from being 
lighter the structural performance offered by it far more versatile than can obtained 
with the conventional materials. Composite plates with holes and cutouts are ex- 
tensively used in aerospace structures. Cutouts change the mechanical behavior 
of coitiposiUi plates. Hence the stability analysis of such plate structures is very 
essential for understanding the behavior of systems under different loadings. 

In the present study an attempt has l)een made to take into account the effects of 
the cutouts in a moderately thick/very thick composite laminated plates under in- 
plane, compressive and shear loading using "Simple Higher Order Shear Deformation 
Theory" based on four unknown displacement functions instead of five, which is 
commonly used for most of the other higher order theories.The finite element method 
is employed to study the buckling load. A C* continuous shear flexible finite element 
based on the HSDT model is developed using the Hermite cubic polynomial for the 
rectangular element. The plate geometry has been modeled using linear Lagrange 
rectangular elements and hence subparametric formulation has been adopted. 

The analytical results obtained for a plate without a cutout have been compared 
with the results obtained by different theories. The change in buckling response of 
thick rectangular laminated plates with and without cutout with respect to the fibre 
orientation angles has been studied. The effect of the cutout dimensions to the plate 
dimensions on the buckling load is plotted for cross-ply laminate. The interaction 
curves for and Ny and for and N^y with and without a cutout have been 
studied . 
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Chapter 1 
Introduction 


Structural materials can be broadly divided into four basic categories:metals, poly- 
mers, ceramics and composites. The word composite literally means consisting of 
two or more distinct parts. Thus a material which consists of two or more distinct 
materials combined in a macroscopic scale is called a composite. Common materi- 
als which can be included in this broad definition include reinforced concrete, fibre 
reinforced plastic, plywood, etc. Among the two phases one acts as as a binder or 
matrix and the other acts as a reinforcement. 

Composites are generally used in structural components because they have de- 
sirable properties which could not be achieved by either of the constituent materials 
acting alone. The most common example is the fibrous composite consisting of 
reinforcing fibres embedded in a binder, or matrix material. Particle or flake rein- 
forcements are also used but they are not so effective as fibres. Thus, in describing 
a composite material as a system, besides specifying the constituent materials and 
their properties, one needs to specify the geometry of the reinforcement with ref- 
erence to the system. The geometry of the reinforcement may be described by the 
size, shape and size distribution. 

Broadly composites can be classified into two categories: 

1. Particle-reinforced composites. 

2. Fibre-reinforced composites. 
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In the present investigation laminated fibrous composite is considered. Some of 
the commonly used fibre materials are E-glass, S-glass, Kevlar , Graphite, Boron, 
Silica, IHingsten etc. h'ibres are us(!d to provirled the necessary strength and stiffness. 
Matrix materials (Polymers) can be classified into two types depending on their 
structure and thermal behaviour. Those are namely: 

1. Thermoplastic Polymers -Soften and melt on heating. Polyethylene, Polystyrene, 
Nylon, Placentals belong to this group. 

2. Thermosetting Polyrners-Do not soften but decompose on heating. Once so- 
lidified by a cross-linking curing process, tlusy cannot be reshaped. Common 
examples of this kind include Epoxies, Polyesters, Phenolic, Silicone etc. 

The main advantage of Polymeric matrix are low cost, easy processibility, good 
chemical resistance and low specific gravity. 

Composites offers unique opportunities in design. Apart from being lighter it is 
a good substitute for conventional metallic materials. The structural performance 
offered by it is far more versatile than can be realized with conventional materials. 
An in-depth understanding of the principles governing their structural behaviour 
can help us explore this versatility to a maximum possible extent. 

1.1 Composite Plates with Cutouts 

Composite plates with holes and cutouts are extensively used in aerospace structures. 
The holes and cutouts are made for saving weight, for accommodating fuel, hydraulic 
or electrical lines and for providing accessibility to other parts of structures. Most 
of the studies in the stability analysis of plates are dealt with plates which are 
continuous or complete. In reality most aerostructures are flawed (either a crack or 
a cutout) in some way. Cutouts change the mechanical behaviour of plates. Hence, 
the analysis of such plate structures is mush essential. Thus the stability analysis of 
structures with cutouts is of technical importance for understanding the behaviour 
of systems under different types of loadings. Despite the practical importance of 
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these structures with cutouts, the number of technical papers and reports dealing 
with the subjects are very limited due to the complexity involved. 

1.2 Motivation 

'i'he key factor in utilizing the strength and uniciueness of laminated composite 
plates is the proper understanding of its structural response under different work 
load conditions. Plate-like subcomponents come in various forms such as a annular 
plate or a rectangular plate that has a cutout. The present study focuses on the 
rectangular or square plates that have cutouts at different locations. Developing 
a thorough understanding of the behaviour of this subcomponent is a fundamental 
step towards understanding the behaviour of complex structures with cutouts such 
as airplane wing ribs, Knowledge of the basic response of the subcomponent provides 
useful information for the preliminary design of complex structures. In addition, this 
basic knowledge provides valuable insight into modelling complex structures with 
general purpose finite element codes, a step that takes place at a later phase in 
design process. 

One behavioral aspect of the laminated composite plates which is vital for safe 
design is the initiation of buckling under the action of in-plane compressive and 
shear loads. Buckling can be termed as an unstable equilibrium of the body. The 
load at which an initially flat configuration of the plate under the action of in-plane 
forces suddenly changes to a bent configuration is known as initial buckling load. 
As the structure undergoes failure at this point it is imperative that a considerable 
arriouiJt of attention liiust be paid to find out its initial buckling load. The exact 
determination of this-'load governs the safe design process. 

Most of the advanced composites in use have low transverse shear modulus (due 
to polymer matrices) compared to their in-plane moduli. Therefore, the classical 
theories of plates which predict the response characteristics of thin plates fairly 
accurately, may not yield accurate results for the similar configurations made up of 
composite materials. Thus, transverse shear deformation plays an important role 
in reducing the flexural stiffness of laminated composite plates. The classical plate 



thoory due to Poisson-Kirchoff neglected the transverse shear flexibility and assumed 
inextensible straight normals.This leads to an under-prediction of deflections and 
over-prediction of buckling loads. 

So, for a better estimation of the initial buckling of composite laminated plates, 
many theories have been proposed over the years which include the effect of trans- 
verse shear deformation. The First order shear deformation theory due to Mindlin 
for thick plates assumes constant transverse shear strains through the thickness. As 
this theory allows only constant shear strain and does not satisfy the stress free 
conditions at the top and bottom surfaces of the plate, appropriate shear correc- 
tion factors are introduced to achieve meaningful results. The exact evaluation of 
these shear correction factors is very difficult as it depends on the geometry and 
Poission’s ratio. In recent years several simplified higher order shear deformation 
theories (HSDT) have been proposed, which not only allow parabolic variation of 
transverse shear strains but also satisfy the vanishing of transverse shear /strains at 
the bounding planes of the plate have been proposed. These higher order theories 
normally reepure five unknown displacement functions and requires continuous 
elements for finite element analysis. This requires numerical strategies such as re- 
duced/selective integration to avoid shear locking effect (see [1]). A higher order 
theory proposed by Lim et al. [2| involves only four unknown displacement functions 
which requires continuous plate bending element and therefore does not neces- 
sitate numerical strategies to avoid shear locking effect. Also this theony involves 
lesser number of unknowns and is expected to be computationally more eflFective 
than the other theories which involve more unknowns. 


1.3 Literature Survey 

The advent of new strong and stiff fibre-reinforced composite materials approxi- 
mately four decades ago revolutionized the structural applications in the aerospace 
industry. Since then, plates made of such materials are being extensively used in 
the industry. As a result various theories have been developed to understand the 
mechanical behaviour of composite laminates. Some of the studies ( [3] and [4]) have 
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shown that the transverse shear ed'ect is quite significant in the layered composite 
plates due to high ratio of in-plane elastic modulus to transverse shear modulus 
which made classical laminate jdate theory unsuitable for the analysis. To over- 
come this difficulty a modified bending theory was formulated which accounts for 
the transverse shear deformation in isotropic plates. This plate bending theory 
is known as Reissner/Mindlin plate bending theory. This theory was applied for 
anisotropic plates by Whitney and Pagano (see [3] and [4]). In this theory the 
transverse shear strains are constant through the thickness of the plate. A shear 
correction factor has to be introduced to account for the discrepancy between the 
theory and the actual behaviour. Though this theory gives better results than clas- 
sical laminate theory, the error in the solution increases as the plate thickness to 
side ratio increases. Thus to account for a better representation of the shear distri- 
bution across the thickness, Reddy [5] proposed a higher order shear deformation 
theory. A simplified higher order theory proposed by Lim et al. [2] involves only 
four unknowns instead of five for Reddy |5]. The simplified theory allows for the 
use of continuous plate bending element which is free from any shear locking 
effect. Singh [6| and Chakraborty [7| used this model for the response of laminated 
composite plates and beams. The stability of plates with perforated plates is stud- 
ied by Yettram et al. [8|. They employed conjugate/load displacement method to 
analyse the structural stability of the perforated plates. They made the analysis on 
isotropic plates. Srivasta and Murti [9j presented a parametric study of the com- 
pression buckling behaviour of stress loaded composite plate with a central circular 
cutout. Jain and Kumar [10] analysed the postbuckling response of square lami- 
nates with a central circular/elliptical cutout. In their study, they investigated the 
large deflection behaviour and postbuckling behaviour of axially compressed square 
laminates with a circular and elliptical cutout. Prabhakara and Datta [11] studied 
the vibration and buckling behaviour of plates with centrally located cutouts. To 
the best of author’s knowledge no results are available for the initial buckling anal- 
ysis of laminated plates with rectangular cutout under uniaxial, biaxial compressive 
and shear in-plane loading in the published literature, using the simple higher order 
deformation theory. 
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1.4 Present Investigation 

In the prcHont investigation tlie following studies have been done. 

1. Buckling behaviour of laminated composite plate without a cutout under in- 
plane compressive loading using simple' liighe.r order shear deformation theory 
as mentioned in section 1.3 is studied. The results are validated against those 
of other investigators. 

2. Buckling behaviour of laminated composite plate with cutout under in-plane 
loading using simple higher order shear deformation theory. 

3. The effect of various parameters like cutout size, plate aspect ratio, length to 
thickness ratio and fibre orientation on the initial buckling load of laminates. 

4. Development of interaction curves for biaxial loading {Nx and Ny and Nx and 
+Nxy, Nx and -Nxy) of symmetric laminates for different cutout sizes, aspect 
ratio, length to thickness ratio and lamination angles. This will help in proper 
choice of laminate dimension. 

1.5 Layout of Thesis 

The thesis is broadly divided into five chapters.Thc details of thesis layout are pre- 
sented here. 

1. The first chapter is devoted to the the motivation for the present work, liter- 
ature review, the focus of the present study and the layout of the thesis. 

2. The second chapter discusses the plate model using higher order shear de- 
formation theory, formulation of the problem and the laminate constitutive 
equations. 

3. The third chapter deals with the finite element formulation of the present 
problem using energy principle. 
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4. Validation of the present model and several numerical results that detail the 
present investigation are presented in the fourth ehai)ter. 

5. The fifth chapter deals with the analysis of the results, concluding remarks 
and the further extension of the problem. 


9 



Chapter 2 

Higher Order Plate Theory 


2.1 Introduction 

The classical theory of plates, in which it is assumed that normals to the mid- 
plane before deformation remain straight and normal to the plane after deformation, 
under-predicts deflections and over-predicts buckling loads. This is because the 
transverse shear strains are neglected in the classical theory. Many plate theories 
exist that account for transverse shear. Of these, the theories based on the assumed 
displacement field provide a background for the ijresent theory. This chapter deals 
with the details of the plate theory proposed by Lim et al. |2] for isotropic plates 
and is extended to laminated composite plates. 

The displacement field includes classical plate theory and first order shear defor- 
mation theory as subsets and accounts for paral)olic variation of the transverse shear 
strains as well as the surface boundary conditions of zero transverse shear stresses 
(and hence shear strains) at the top and bottom surfaces of the plate. Hence there is 
no need to use shear correction coefficients in computations. The displacement field 
proposed by Reddy [5] is modified by reducing the number of variables to four (five, 
in case of Reddy’s [5] theory). This simplification is based on the assumption that 
in-plarie rotation tensor is constant through the thickness and the transverse dis- 
placement can be separated into two components in such a way that the transverse 
shear strains are functions of only one of these components. 
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Figure 2.1: Laminated plate with orientation of global co-ordinates 


2.2 Formulation of Displacement Field 

The displacement field shown in Fig. 2.1 u{x, y, z), v{x, y, z) and w{x, y, z) is repre- 
sented as: 


u{x, y, z) = iLo{x, y) - zwb,^[x, y) + z‘^(().j,{x, y) + z^'tpx{x, y)\ 

v{x,y,z) = VQ{x,y) - z%Vb,y{x,y) + z‘^(t)y{x,y) z^'iljy{x,y)\ ( 2 . 1 ) 

w{x, y, z) = wo{x, y) + Ws{x, y) 

where uq, vq and w denote the displacements of a point (x, y) on the mid- 
plane, and 'tpx s-nd 'ipy are the warping of the normals to the midplane about the Y 
and X axes, respectively.The transverse displacement component Wb is such that its 
derivatives are numerically equal to the rotation of the cross-section (i.e. (j) = — viUi,) 
and Ws is the displacement due to the effect of transverse shear deformation of the 
cross-section. The assumjjtion that Wb and are functions of x and y only is 
justified, since the transverse normal stress is of the order of {h/a)'^ times the in- 
plane normal stresses. 

The von-Karman type nonlinear strain-displacement relationships for the higher 
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order theory under consideration can bo written as follows: 


^xx d" d* d” 

^VV "h 2 {^^biy d" “ Z'^^btyy d" z 4^yiy d" z' VV’yi 

e^z = 0; 

"7j:y ^k)yy I d' 1' ti lx ) ( I lUfljy) 2z'tll/i,.i^y (2'2) 

d''2^' {.4'yix d" 4^Xyy) d" z {‘l/jy,^ + i^Xty)] 

Ixz = Wsrx d- 220X d- ‘iz^i>x\ 
lyz = Ws,y + 2z(l)y + 


whore coinina (,) denotes the partial derivatives. 



Figure 2.2: Representation of arbitrary three dimensional domain 


The surface boundary condition that the transverse shear stresses vanish on the 
plate’s top and bottom faces (Fig. 2.2) is equivalent to the requirement that the 
corresponding strains be zero on these surfaces, i.e. 

'Yxzi^i Vt d-^) ~ Ut 0 

On introduction of the conditions as given above in the expressions for transverse 
shear strains ( from Eq. 2.2 ), the following relations are obtained 

4>x ~ 4’y ~ ^ and 
i^x — ~^'bb>a>x ) V'v ~ 
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I'tu! (lisplaccnient field of E(}. (2.1) is modified by setting (p.j. and cf)y to be zero 
according to conditions of E(]. (2.3). The resulting displacement field is given as : 

u{x, y, z) = uoix, y) - zwb,^(x, y) - |^w;«,^.(.x, y)\ 

v{x, y, z) = vo{x, y) - zwi„y{x, y) - y)\ (2.4) 

w{x, y, z) = wo{x, y) + v)^{x, y); 

In Eq. (2.4) n, v and w arc the displacements along x, y and z directions re- 
spectively. uq,vq and w (twi, and u;.,) are the mid-plane displacements. Thus, the 
generalized displacement vector {d} of the mid-surface contains four degrees of free- 
dom (DOF) and is given by; 


{ 5 } = {uo,Vo, 

The corresponding strain-displacement relationships are: 

£xx = '«0,x + - 5P'“'.,,xxi 

hV = - ZW,„yy - 

Czz = 0 ; 

7x2/ = UQ,y + + w^,.J{wi„y -I- ■w,„y) - 2zxvb,^y (2.5) 

— ^w 

7xz ~ (f TT^ )^^'’a)X) 

7j/ 2 = (1 “ ■jpr)'^sij/i 

It can be noted that these kinematic relationships can be obtained by substituting 
■0X = -Wbyj;, xpy — —Wb,y and w = Wb+Wa in Reddy’s [5] higher order theory involving 
five unknowns (uf), Uo, ujo, r/^xi V'y)- 

It can be shown that the present formulation degenerates to first order shear 
deformation theory by suitably dropping higher order terms and incorporating ap- 
propriate shear correction factors (for isotropic plate whose value is 5/6). The 
classical laminate theory (CLT) can be derived from the present theory by equating 
Wa to zero. 
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The strain-clisplacenieut relationships given in Eq. ( 2.2) can be rewritten as: 


£‘.rx — COx + 

^yy — ^Oy + ZK'Oy + Z^Kly] 

^zz = 0 ; ( 2 . 6 ) 

Ixy = TOxy + ZKoj;y + Z^/t/x-yi 

Ixz = (1 

Az^ 

Tyz (1 y^2 

wliere : 

^Ox “• ^Oy ^ ‘^V)jy! 70x-y 

^Ox’ ^^b)xx) ^Oy ^^b)yy) ^{)xy 


^0)x) 

^ 

3^2 


2.3 Laminate Constitutive Equations 

A single layer of a laminated composite material is generally referred to as a ply 
or lamina. A lamina (considered a unidirectional composite) is characterized by 
having all fibres (either a single ply or multiple plies) oriented in the same direction. 
This model allows one to treat the lamina as an orthotropic material whose material 
symmetry planes are parallel and transverse to the fibre direction. The model used 
to represent a lamina consists of unidirectional fibres per layer which is further 
assumed as perfectly straight and uniformly oriented within the lamina. 

The following figure represents a laminate with the material co-ordinates axes L 
and T defined parallel and perpendicular to the fiber direction respectively. X, Y and 
Z represents the global co-ordinates. The material co-ordinate axis T' is assumed 
to be along Z axis. The angle between global co-ordinate axis X and material co- 
ordinate axis L is 6 and is known as the fibre orientation angle. As a sign convention 
anti-clockwise angle from -j-ve X direction is taken as positive. 
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Figure 2.3: Co-ordinate axes for the lamina 

2.4 Generalised Hooke’s law 

In the derivation of the laminate constitutive eq\iations the following two assump- 
tions are made. 

1. The lamina is a continuum. 

2. It behaves as a linearly elastic material 

Further at the micro-level the following assumptions are made about the mate- 
rial. 

1. Perfect bonding exists between the fibres and matrix. 

2. Fibres are parallel and distributed uniformly throughout the laminate. 

3. The matrix is free of voids or micro cracks and is initially in a stress free state. 

4. Both fibres and matrix are isotropic and obey Hooke’s law. 

Stress-strain relations for the lamina in the material co-ordinate axes, whose 
fibres are oriented at an angle 6 with reference to the X axis is given as (see[ll]): 

{cnh = [QijUcj}k (2.7) 


15 




wlun-e {(Til iH the vector of stress components, [Q,j]fc is the 2D- stiffness ma- 
trix with respect to material co-ordinate axes, and {cj} are the engineering strain 
components for the lamina. 

In the expanded form, the above relation can be written as : 


/ \ 
ctl 


(JT 


Tir 

> = 



\ ^ ✓ 

k 


Qn 

Qi2 

0 

0 

0 


( ^ 

Q 21 

Q 22 

0 

0 

0 


ct 

0 

0 

Qm 0 

0 

< 

ILT / 

0 

0 

0 

Q 44 

0 



0 

0 

0 

() 

Q(j(j 

k 

, 7rT' , 


where ; 


Qu = 
Q22 = 


(1 - 

Et 

(1 - ’ 

Q44 = ^ lt' > 


Qi2 = Q 2 I = 


Ex^it 

(1 - UiT^^TlV 


= Glt\ 
Q55 = G ^' x '] 


(2.8) 


The stress and strains in the X, Y and Z directions are obtained by the trans- 
formation of the relations given in Eq. (2.7). The transformed stress strain relations 
for the kth lamina are given as; 


/ \ 

^ X 


Qn Qi2 Qig 0 0 


r \ 



Qi2 Q22 Q2G 0 0 


Gy 

< Txy 

> = 

Qi6 Q2G Qgg 0 0 

< 

'Ixy ^ 

'^XZ 


0 0 0 (344 Q 45 


'Ixz 


k 

0 0 0 (345 Q55 

k 



(2.9) 
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\vh(‘r(i: 


On = Oi 1 0 + ( 2 Qn + + Q22A 

Qi2 = Qi 2(^‘^C^) ■+ (Qn + Qi 2 - 4QG6)il'^C^; 

O16 (On ~ Qi 2 + 2Q{i(i)6’c^ + (O12 ~ O22 + 2Q(j0)s^c; 
O 22 = On*'"* + ( 2 O 12 + 40(iG)4'''^(-^ ^■ 022^-^ 

O26 = (On - O12 - 20 cg)s^c+ (O12 - O22 + 20og)sc^; 
OoG = (On ~ 2O12 + O22 — 20 gg)s^<^^ + Ogg(c^ + s'*); 

044 = 044 C^ + OssS^i 

045 (O44 Oss)^®) 

O55 = Oas^^ + 044 S^; 


where, [Oij] is the transformed stiffness matrix, c — cosff and s = sin 0 with 0 
denoting the angle of orientation of the kn^ lamina. 

Using the above lamina constitutive equations and integrating the stresses over 
the laminate thickness, the stress resultants, stress couples, transverse shear resul- 
tants and higher order stress couples per unit length, for a laminate with NL laminae 
are, 


' {Ni} ' 


M«1 (S«l [E«l |0| 



{Mi} 

^ = 

Wil ! 0«1 [P«l [o| 

< 

{kOj} 

i 

{Pi} 


lEij] (Rj) m ioi 


{%} 

{Q.} J 


. (01 (0| |(>1 |(7„„1 . 


^ {77/1} , 


(i , j = 1 , 2 and 6 ) 
(1 , m =: 4 , 5) 


( 2 . 10 ) 


(Detailed expansion of Eq. (2.9) is given in Appendix A.) 

where {gq} are the mid-surface strains, {kq} are the mid-surface bending and twisting 
curvatures, {k;} are the higher order terms and { 7 } are the transverse shear slope 
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and 





1 


{Ni}'" 

= (/Vx.,/v„yv,,, 

! 

^ 

Elis 

11 

dz-, 

1 

(2.11) 

{Mip = 


’=sy 

*,■=1 

ip^xt ^y-i '^xy) ^ 

(2.12) 


= (C,, Py, P,y) 

k=i 

( ^x ) ^y ) T'xy ) 

(2.13) 


iV Li I y1 2 

{Ql}’’ = IQ,., ft) = E / ' (1 - %) dz; (2.14) 

fc=l 


Here, Zi are the x coordinates corresponding to the laminae interfaces as shown in 
Fig. 2.4. 



Figure 2.4: Geometry of multi layered laminate 

Thus, with the assumed displacement model, the various stiffness matrices derived 
are: 

[Aij] = Extensional; [By] = Bending-extensional coupling; 

[Dy] = Bending; [G;r„] = Transverse shear stiffness and 
[Ey], [Fjj], [Hij] ~ Higher order stiffnesses. 

The set of stiffness matrices [Aij], [By], [By], [Ey], [Fy], [By] and [Gj^] are used in 
forming overall stiffness matrix [B^] for the laminate (for details see [7]). 
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Chapter 3 

Finite Element Formulation 


This chapter deals with the finite element formulation of the problem under consid- 
eration. Four noded rectangular elements are used to represent the total domain of 
the plate. To satisfy the continuity of the slope, Hermite cubic interpolations are 
used in the approximation. A continuous shear flexible element based on the 
presented higher order theory is developed using the Hermite interpolation formulae 
as indicated by Bogner et al. [13]. The geometry of the plate is modelled using the 
linear Lagrange interpolation functions. A typical mesh generated over the plate 
domain with and without a centrally located cutout is shown in Figs 3.1. and 3.2. 
The details of the finite element formulation is presented below. 


3.1 Definition 


Let V be the displacement vector defined as: 



(3.1) 


The stress and strain vectors corresponding to {F} be {cr} and {e} and can be 
defined as; 
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Figure 3.1: Typical mesh of Rectangular plate domain 



From the generalized Hooke’s law which relates the stress components to the 
respective strain components in global coordinate system: 

{a) = (f51 (£) (3-3) 

The material stiffness matrix [Q] tor ortliotropic material is given in Eq. (2.8) tor 
each lamina. 

We are concerned with the initial buckling of the laminate. Hence, the non-linear 
terms in the strain-displacement relationships as described in Eq. (2.5) are neglected 
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c is the length of the cutout 
d is the width of the cutout 

Figure 3.2; Typical mesh of Rectangular plate domain with a cutout 


for our study bcicause the transverse d(!flection (both ujf, and w,, ) is assiniicd to be3 
very small. 

The linear strain-displacement relationships can be written as; 


^xx — ^ 

^yy ~ ^ 

= 0 ; 


4*3 

57 ? 





8 

Jxy = '^^Oiy + “ 2 Z Wyyj,y — 

ixz= (1 “ 

7v*= (1- 


The above relations are expressed in matrix form as ; 


(3.4) 


{e} = [D]{V} 


(3.5) 
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where [D] is the differential operator in terms of global coordinates, defined as 




d 

dx 

0 

0 

0 

Jl 

dy 

0 


Ji 

0 

dy 

dx 

d ' 
Dz 

0 

9. 

dx 

0 

d 


dz 

dy 


(3.6) 


The components of displacement can be expressed in terms of the four unknowns, 
which can be written as; 


{Sy' = {u„,vn,Wb,Ws} 


(3.7) 


For any element the field variables can be written in terms of shape functions 
and nodal variables as: 

4 4 ; 4 

^0(2;, y) = ^ '*'■<), xi + X/ 


1=1 

4 


t=l 


i=l 


i=r 


i=l 

4 4 

~ ^ ^ -^1+4 ^Ojari ^ ^ 

t=l i=l 

4 

^1+12 
i=l 

4 4 

1=1 1=1 

4 

^ V N^4.12 '^h)vxyi^ 
i=l 

4 4 

Ws{x,y) = X + X' S ^^+8 

1=1 

E"« 


(3.8) 


1=1 


i=l 


i=l 


'^i+12 


i=l 


In the above expression uoi, ua,xi^ uo,y^,tio,xyii ^Oii 2.'o,y^) ■wonyi) 

refc,^,, tn.i, w,,^y. {i = 1,4) are the sixteen degrees of freedom 
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per node, iV-s are the element shape functions in local co-ordinates (in terms of 
X and Y as shown in Fig.3.2). In natural co-ordinates (in terms of ^ and rj as 
shown in Fig.3.2) iV(s can be expressed as (see [1]). Substituting Eq. (3.7) in 



Four noded rectangular element in Four noded square clement in 

local co-ordinate system natural co-ordinate system 


Figure 3.3: Four noded rectangular element 

the strain-displacement relations given in Eq. (3.4) the following relations can be 
established: 

U)x 

70 ;rj; 

^Ox 

^Ix 
f^ly 
f^lxy 
'^s ^ 

Wo 

where: 
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{(5} { ^Oi) ^0)xj) ^Oixyii ^Oixy^j ^fcji '^b)xi^ '^bixj^i 

and: 

^O.T ^'Oix) f()y ~ ^0’j/> TOxy ~ ^Oiy d" ^Oix! 

^Ox ^6)xx) "^^djj'yy) ^'Oxy 

^Ix ~~ I^ly — ~3X5 '^S'yyi ^Ixy — ~3^ '^siiy 

The details of [jBq] are given in Appendix B. 

3.2 Finite Element Formulation 

An eigen value buckling analysis is a linearized form of a geometrically nonlinear 
formulation, useful for estimating the limits of elastic stability. 

The undeforined configuration of the plate is denoted by Q and its boundary by 
dCt. The infinitesimal strain is defined in terms of the Cartesian components of the 
displacements {ui,i = 1,2,3): 

1 . 

which is a simplification of the Green-Lagrangc strains defined by 

“ij — ^ij 2 ^ j 

The simplification is justified by the assumption that << 1 and hence the 
product terms Ua,iUa,j are negligible in relation to The stress-strain relationship 
is: 

O ' ij “ “b C^ijkl ^kl 

where a°j is a pre-existing stress state, independent of Ui, and Cijki is the tensor 
of the elastic moduli of the material. An important property of afj is that, it is in 
equilibrium with the corresponding tractions 7}^ = cr?- rij in the sense: [see [14]j 

i f aUvuj + vj,,)dV = / If Vi dA for all Vi e E {Q) 

^ Jn Jon 
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where dV and dA represent the difl'crential volume and diflerential area respectively, 
and E{Q,) is the space of kinematically admissible perturbations. When the reference 
configuration is stress-free (i.c., crf^ = 0) then the potential energy is defined by: 

(■'0 - ^ f Uj Cfci dV - f Ti Ui dA (3.10) 

^ Jn Jan 

The exact solution minimizes 11 on the sot of all kinematically admissible functions 
denoted by E (Cl). When the reference configuration is not stress-free thou the work 
done by a?- due to the nonlinear strain terms may not be negligible. Therefore the 
potential energy expression is written in the following form. 

n(u) = \ I Cijki eij Cki dV + ^ / dV - [ TiUidA (3.11) 

2 Jn 2 Jii Jqq^ 

The first integral in the above eijuation represents the internal strain energy and the 

second term represents the work done by the initial stresses due to the nonlinear 

strain terms. The work done by cr?. due to the linear strain terms is canceled by the 

work done by Tf in the sense of Eq.(3.10). 

The exact solution to this problem is obtained by minimizing the total potential 

energy 11. This can be obtained as: 

II = 0 (3.12) 

which is also the Viiiual Work Fonnulation of the problem in terms of the com- 
ponents of ^(and hence 6). This leads to the generalized finite element formulation 
which can be expressed as: 

1K|(5} = {F) (3.13) 

Where [i?] is the global stiffness matrix which is composed of two components (i.e. 
global elastic stiffness matrix [K] and global geometric stiffness matrix [/f^]) and 
{F} is the global load vector which is zero in our case since = 0. 

Hence Eq. (3.15) can be rewritten as: 

’[/^1 - UQ ]{' 5 }= 0 ( 3 - 14 ) 

In the following sections, we are going to employ the virtual work formulation 
to derive the finite element equations governing the initial buckling problem. 
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3.3 Deriving the Element Elastic Stiffness Matrix 


The elastic stiffness matrix corresponding to .assumed do.fonnation state of an ele- 
ment can be defined by expressing the internal strain energy in terms of unknown 
nodal displacements. In the formulation of unsyrametric laminates the membrane, 
the flexure, membrane-flexure coupling and shear strains contribute to strain energy. 
The internal strain energy of an element (given by area A^) can be expressed as: 

O / ^ ^ 0 . 7 : 3 / d“ Alx T T AIxy I^Oxy d” ^Ix 

T ^ly d" ^ xy ^Ixy d" T ^'A (3.15) 

The Eq. (3.8) along with Eq. (2.9) can be used to express the internal strain 
energy as; 

t'" = 5 / ( {i'F IBo'l [0.1 [Sol (i'}) dA (3.16) 

The strain energy expression can be written in concise form as: 

lyM = 1 {5^} (3.17) 

where [K^] is the elastic stiffness matrix and {5'’} are the unknown displacements 
corresponding to the finite element solution in an element e . [K’’-] is given by: 

[BoF [D.I [Sol dA (3.18) 

Numerical integration is carried out to get the element elastic stiffness matrix. 

3.4 Deriving the Element Geometric Stiffness Ma- 
trix 

The total work done by the in-plane applied loads in an element can be expressed 
as: 

T' = 5 / { {«“}’■ [Bjl m 1B,1 {«'}) dA (3.19) 

T® can be written in a concise form as; 

T{«) = 1 f (3.20) 
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The details of [/i^] is given in Appendix C. 
whereas: 


[N] = 


N. 




N N 

J-yxy iyy 


(3.21) 


Assuming that forces Nx, Ny and Nxy are represented by certain expressions with 
a common factor A, so that 


then: 


or 


Nx -XNx, N„ 

[A] 

[A] = A 


A Ny, - A Nxy 


Nx 

Nxy 


Ny 

(JV| 



(3.22) 


where: 


[N] 


Nx Nxy 
Nxy Ny 


(3.23) 


where [A^] is the geometric stiffness matrix and {(f®} are the unknown displace- 
ments corresponding to the finite element solution in an element e. 

[A®] is given by: 

[I<1] = A / [Byf [A] [B,] dA- (3.24) 

or 

(/fg = A [f?/| (3.25) 

Numerical integration is carried out to get the clement geometric stiffness matrix. 

3.5 Methodology for Finding Initial Buckling Load 

The main objective of the present analysis is to investigate the initial buckling load 
of the composite laminate under different boundary conditions, cutout sizes, location 
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of the cutout, fibre orientation angles, plate aspect ratio’s etc. Since we are dealing 
with the physical stability of the system, this buckling problem is an eigen value 
problem which is described below. 

From Eqns. (3. 16), (3.20) and (3.23) we obtain the final FEM equation as follows: 

■ [K\ - j {J} = 0 (3.26) 

The above equation can be stated as "Assuming that the steady state solution 
of the problem is known, is there any other solution into which the system would 
bifurcate if it were slightly disturbed from it equlibrium position". 

The characteristic equation of the above eigenvalue problem can be expressed as: 


[K] - X[I<,] 


= 0 


(3.27) 


Eq. (3.25) is solved using NAG ROUTINES to obtain the eigenvalues Xi (i = 1,M) 
where M is the dimension of the matrices [K] and [Kg] (and is equal to the total 
number of degrees of freedom of the problem). The minimum value of A corresponds 
to the scaled initial buckling load from which we obtain the initial buckling load as 
follows: 

Kx ( or Ky or Kxy ) “ ^min Njp ( or Ky or Kxy ) 

The corresponding eigenvector (which is also found using NAG ROUTINES) gives 
the mode shape (not the exact displacements) corresponding to the initial buckling 
load. 
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3.6 Geometric Approximation 

3.6.1 Linear Mapping of Straight Edge Rectangular Elements 

In any finite element formulation an accurate representation of the domains is very 
essential. Any domain can be represented using rectangular elements since the inter- 
polation functions are easily derivable for them. Numerical integration schemes such 
as the Gauss-Legendre scheme require that the integral he evaluated on a specific 
domain or with respect to a specific co-ordinate system.The transformation of the 
geometry and the variable coefficients of the differential equation from the problem 
co-ordinates ( x, y) to the co-ordinates (^, ?;) results in algebraically complex expres- 
sions and this precludes analytical (i.e., exact) evaluation of the integrals. Thus, the 
transformation of the given integral expressions, defined over the element f2 to one 
on the domain D, must be such as to facilitate numerical integration. Each element 
of the finite element mesh is transformed to Tl only for the purpose of numerically 
evaluating the integrals. The element il is called the master element. The geometric 
transformation from the physical element to master element can be done by linear 
mapping which is discussed below. 


Y 


(X4,Y4) (X3,Y3) 



(X2,Y2) 



X 


.- 1 ) 


( 0 , 0 ) 


(- 1 ,- 1 ) ( 


, 1 ) 


.- 1 ) 


Physical Element Master Element 

Figure 3.4: Linear mapping in two dimensional domain 

The transformation between and (l (or, equivalently between {x and y) and 


2 !) 



and 7]) is accomplished by a coordinate transformation of the form 


4 4 

= = 0-28) 
J=1 j=l 

where Xf' and Yj^ are the global coordinates on the physical element nodes and 
ifjj denote the finite element interpolation functions of the master element Cl. Sub 
parametric formulation is used to represent geometry of the domain. Hence bilinear 
Lagrange interpolation functions are employed here, 
are defined in the following way. 

^ (1 - 0 (1 - h ); 

V'a" = I (1 + 0 (1 - h); 

VV = ^ (1 + 0 (1 + v)\ (3.29) 

V'/ = ^ (1 - 0 (1 + r,); 

The transformation can be expressed in terms of the local coordinates ^ and 
t] by means of the. Hence by the chain rule of partial differentiation, we have 

dijji^ _ ' dx dij)i ' dy ^ 

dx ()( dy dl' wL 

^ ^ ^ ^ 
di] dx dr] dy dy 

which gives the relation between the derivatives of ipi^ with respect to the global 
and local coordinates. 

The matrix in the above equation is called the Jacobian matrix of the transfor- 
mation. 

I Transforming dxdy into d^drj 

Let the differential area dxdy is formed through vectors dx and dy with magnitude 
dA and direction normal to the elemental area is K 
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dA = dx dy = [dx y. dy] .k 


(3.31) 


dxdy = [ ( ^ i + ( ^ ) X {^d^i + 


Ot] 


<iv J ) ] 


dx Ox 
drj 


d^ 


dy 


] d^ dy 


= \J\didy 


(3.32) 


where 


.^._dx dy dy dx 
d( dy dy 

Hence f {x, y) dx d.y can be written as 



/ 

crn 


f V)d$dy 


(3.33) 


(3.34) 


such that / ( C, 7?) = \J\ f {x ( ?;), y{ e, y ) ) 

The condition which must be satisfied for physically meaningful mapping of the 
domain is that,| J| > 0. 


I Numerical Integration 

Quadrature formula for integrals defined over the rectangular master element CIr 
can be derived as follows. 


„ M N . 

/ /(?,>))<(? EE / /(«., Wi Wj (3.35) 

Jd i.i M 

where M and N denote the number of quadrature pints in the ^ and y directions, 
(4ii yj) denote the Gauss points, and Wi and Wj denote the corresponding Gauss 
weights. Order of the' numerical integration should be greater than or equal to (p + 
l)/2 where, p is the order of the integrand. 

In the present formulation 5x5 Gauss quadrature points have been used. 
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Chapter 4 


Numerical Results and Discussion 


4.1 Introduction 

The present chapter deals with the validation and the various results obtained using 
the finite element formulation described in the previous chapter. The main aim 
is to study the initial buckling behaviour of laminated composite plates with and 
without cutouts, using a simple higher order shear deformation theory. For the 
present study it is assumed that there is no coupling between the in-plane and out- 
of-plane displacements, ie., symmetric lay-ups assumed. If this is not true, then 
even for an infinitesimally small in-planc compressive force the laminated composite 
plate will undergo flexure. For symmetric laminates this condition is automatically 
satisfied but for antisymmetric laminates there is always coupling between in-plane 
and out-of-plane quantities. In the case of antisymmetric laminates we have assumed 
that this coupling eflfect is negligibly small and was ignored to find the buckling load. 
This assumption although seems strong it can be shown that if the number of plies 
is increased the coupling effect dies out rapidly and the solution approaches the 
orthotropic solution (see [15]). 

In the present study the following assumptions are made. 

• The thickness of each layer is constant. 

• There is no delamination/slip between the layers. 
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• No micro-buckling takes place before the initial buckling. 

• There is no fibre breakage or matrix cracking before the onset of initial buck- 
ling. 

In the current investigation the following nomenclatures are assumed to denote 
the thickness of laminate. 

• ajh > 100 : Thin laminate . 

• 20 < a/h < 50 : Moderately thick laminate . 

• a//i < 10 : Very thick laminate . 

This classification is tentative and actual definition depends on a variety of fac- 
tors such as loadings, boundary conditions etc. 


4.2 Boundary Conditions 


The following types of geometric boundary conditions are considered in the numer- 
ical analysis. 


• u 7 ^ 0 and u = 0 along loaded edges ; 


• « = 0 and u 7 ^ 0 along unloaded edges ; 


• Simply-Supported (S) : v)b = W/, — 0 ; 

• Clamped (C) : Wb = Ws = Wb^^ = v>s,n = 0 ! 

where, n = x or y depending on the side of the plate. 


To define the complete set of boundary conditions on all four sides of the plate 
the following scheme has been adopted which is explained through an example, 
sees denotes the following set of boundary contitions: 

Wb = = 0 atx = 0'iy 


Wb 

= Ws 

ii 

.o 

II 

= Ws,y = 

Wb 

= Ws 

= 'Wb,^ = 

= = 

Wb 

= Ws 

= 0 



fli y = 0 V X 

at X = a'^y 
at y = by X 
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4.3 Material Properties 

Computations were carried out for the following material properties. 

I 

Table 4.1: Material properties 


Material code 

EUGPa) 

Mu 

Et 

Glt 

Et 

Et 

Et 


Ml [7]* 

130 

13 

0.5 

0.5 

0.2 

0.35 

M2 [12]* 

200 

40 

0.6 

0.6 

0.5 

0.25 

M3 

100 

isotropic 

0.25 

M3 [7|* 

153 

15.904 

0.616 

0.616 

0.336 

0.32 


^Numbers in the brackets represents the references given in the bibilograply. 

4.4 Convergence Study 

In this section, the convergence study of the laminated composite plates without 
cutouts is discussed. 

Convergence of the initial buckling load with mesh refinement is as shown in Ta- 
ble 4.2. Here for all cases Material Code = Ml, Aspect Ratio = 1.0 and length/thickness 
ratio = 10. 

From Table 4.2 it can be concluded that for a square plate, a 10 x 10 mesh gives 
sufficiently accurate results ( for both simply supported and clamped boundary 
conditions ) from engineering point of view. 

For rectangular plates, it has been observed that for a/b = 2 sufficiently con- 
verged result can be obtained with a 10 x 5 mesh whereas for a/b = 3 a 15 x 5, mesh 
is sufficient for both simply supported and clamped boundary conditions. 
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Table 4.2: Convergence study 


Boundary 

condition 

Laminate lay-up 
sequence 

Mesh size 

Nondimensionalized buckling load, 

{N^ = NjyEth^) 

SSSS 

[90/0], 

2x2 

7.4943 

4x4 

9.3565 

8x8 

9.4372 

10 X 10 

9.4441 

CCCC 

[90/0], 

2x2 

12.6457 

4x4 

15.3898 

8x8 

15.9465 

10 X 10 

16.0501 

SSSS 

[90/0/90/0]2 

2x2 

10.3057 

4x4 

10.9218 

8x8 

11.0400 

10 X 10 

11.0477 

sees 

[90/0/90/0]2 

2x2 

24.2386 

4x4 

27.1568 

8x8 

27.9157 

10 X 10 

27.9212 


4.5 Validation Study 

The validation of the finite element plate model is accomplished by comparing the 
results of the following set of problems with results published in the literature. Here 
for all cases material code = M2, aspect ratio - 1.0, V)oundary condition ~ SSSS and 
laminate lay-up sequence = [O/OOjs. The validation study is made for the problems 
defined in Table 4.3. 

The comparative study of the analytical results from the present investigation with 
the results available in the published literature is presented in Table 4.4. 
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Table 4.3: Problem definitions for validation study 


Problem 

number 

Material 

code 

Aspect 

(o/6) 

Thickness 

ratio(a//?.) 

Boundary 

conditions 

laminae lay-up 
sequence 

■am 


1.0 

100 

SSSS 

isotropic 

, P14 

M3 

1.0 

100 

CCCC 

isotropic 

P15 

M2 

1.0 

10 

SSSS 


[0/90]5 

P16 

M2 

1.0 

10 

sscs 


0/90]5 

P17 

M2 

1.0 

10 

cscs 


0/90]5 

P18 

M2 

1.0 

5 

SSSS 


0/90]5 

P19 

M2 

1.0 

5 

sscs 


0/90]5 

P20 

M2 

1.0 

5 • 

CSCS 


[0/90]5 


From Table 4.4 can be observed that the difFeroncc in the results obtained from 
the present study and those of Ref.[12] is less than 2.5% in all the cases. Hence 
we can conclude that the present formulation gives reasonable results and has been 
validated. 
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Table 4.4: Validation study 


Thickness 
ratio ( a/ /i) 

Laminate theory 
[ref. no.] 

Reference 

Present 

Nx 

Ny 

Nx 

Ny 

P13 

CLT [16] 

4.0 

- 

3.997 

- 

HSDT 

- 

- 

- 

- 

P14 

CLT [16] 

10.15 

- 

10.032 

- 

HSDT 


- 

- 


PIO 


35.232 



^36.198 

HSDT [12 

25.828 


25.943 

25.943 

PIG 


- 

59.288^ 

- 

60.563” 

HSDT [12] 

- 

Kkiaatifl 


33.969 

PIT 

CLT [12] 

- 

^^1 

- 



- 

36.657'^ 

- 

iiliiiM:! 

P5 

CLT [12] 

- 

35.232 

- 

36.198 

HSDT [12] 


12.224 

- 

12.603 

■i 

CLT [12] 

- 

59.288'^ 

- 



'mm 

wmm 

- 

12.969 

P20 

CLT [12] 

;Bi 

KBTOiB 

- 

91.015” 

wmmmm 

- 

13.659” 

- 

13.426 


“ Results obtained with the exact solution developed in the cor- 
responding reference. 

’’ Results using the finite elcinent solution in the corresponding 
reference. 

° Results obtained by supressing Wg over the entire domain of the 
plate. 

Nx/Ny are the dimensionless initial buckling values where the 
nondimensionalizing factor is as given in the corresponding ref- 
erences. 
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4.6 Prebuckling Analysis 


The buckling analysis of plates with cutout necessiates the need for inplane stress 
distribution in the plate as it is no longer uniform which is not the case with a plate 
which is continous. Due to possibility of steep stress gradients being present the use 
of a variable mesh spacing is necessary, as it is well known that the success of any 
method for the determination of critical loads is dependent on the accuracy of the 
plane stress distrubution being used. 




Figure 4.2; Graded mesh 
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The in-plane stress resultants are calculated in conventional manner from the 
original displacement functions used in deriving the stiffness of the element. The 
in-plane stress resultants are calculated at the centre of each finite element and these 
are further used in the formation of the geometric stiffness matrix. Two different 
meshes, the equal and the graded mesh, are employed to study the in-plane stress 
distrubution and also the flexure analysis. Figures 4.1 and 4.2 shows the equal and 
graded meshes. The mesh around the cutout is refined in a graded mesh to capture 
the stress gradients near the cutout corners. Figures 4.3 and 4.4 clearly shows the 
in-plane stress distribution in a plate at a given section x-x. From the figures it has 
been observed that the graded mesh gives good results compared to equal mesh. 



Finite elements along Y axis at section x-x 

Figure 4.3: In-plane stress vs Correspondng fi.nite element in Y direction 
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Finite elements along Y axis at section x-x 

Figure 4.4: In-plane stress vs Correspondng finite element in Y direction 

4.7 Validation of Plate Model With Cutout 

The results for an isotropic laminate with cutout is compared with the results given 
in [8]. Figure 4.5 shows the plot of buckling load for an isotropic laminate with 
varying hole ratio with all edges simply supported. From the plot it can be seen that 
the present method gives reasonably good result.Thus a 9 x 9 mesh gives sufficiently 
accurate result for the plate with cutout. It is observed that the buckling load 
decreases continuously with the increase in hole size. For both the finite element 
models the results are very close for small hole size. However the discrepancy is 
observed to be more for larger hole sizes. It is difficult to comment on the finite 
element model in the absence of closed form solution. 


40 





Figure 4.5: Buckling load vs Plate to hole ratio Y-direction under in-plane 
load Ny , (Isotropic Material), a/b=l 

4.8 Parametric Study of Composite Laminates With 
and Without Cutout 

In this section the effects of plate aspect ratio (a/6), length /thickness ratio {a/h) 
and fibre orientation angle {6) on the initial buckling load of laminated ( 8 layer 
antisymmetric, [±^]4 composite plate, with and without a cutout, is studied (using 
HSDT) keeping the material properties constant (Material - M2) under uniaxial 
in-plane compressive loadings. The geometric boundary conditions considered are 
simply supported on all sides. 

The results are presented for the composite laminates with and without a cutout. 
It can be seen from the Figs. 4. 6-4.9 that for a thin laminate (a/h = 100) the fibre 
orientation angle for the maximum value of the fundamental initial buckling load 
( henceforth the fibre orientation angle for the maximum value of the fundamental 
initial buckling load is denoted as {Ocrir.) and the corresponding initial buckling load 
as the critical buckling load) is around 45° for aspect ratio’s (a/6 = 1,2,3). The 
cricical buckling load decreases very slightly with the increase in aspect ratio. 

With increase in the thickness of the laminate transverse shear plays a dominant 
role in the buckling behaviour of the laminate. For a/h = 50 the fibre orientation 
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angles for critical initial buckling load are about 45°, 42° and 40° respectively, for 
aspect ratio’s a/h = 1, 2, 3. 

As the laminate becomes thicker there is a shift in the cricital fiber orientation 
angle towards lesser values. Also this shift is more as the aspect ratio increases. For 
thick laminate {a/h = 10) the effect of fibre orientation angle on the value of critical 
buckling load decreases with the increase in aspect ratio. The decrease in the values 
of initial buckling load is more pronounced for thick laminates. This brings out the 
fact that for thick laminates one should consider transverse shear effects. 

The results presented in Figs. 4.10 - 4.12 show the effect of a centrally located 
cutout on the behaviour of critical buckling load. It can bo seen that the critical 
fibre orientation angle (0cric.) is still 45° for a thin laminate {a/h = 100) for aspect 
ratio {a/h = 1). But for aspect ratio’s [a/h — 2,3) the critical fibre orientation angle 
is 52° and 55° and there is a shift of 10 degrecis, as a/h ratio increases from 1 to 3. 

Figs. 4.12-4.20 shows the variation of buckling load with and without cutout for 
aspect ratio’s {a/h = 1,2,3). From Figs. 4.12-4.20 we could observe that for aspect 
ratio a/h = 1, the critical buckling load of the plate with a cutout approaches that 
of a plate without a cutout for orientation angles C5° to 90°. But this trend is not 
observed for the aspect ratio’s 2 and 3. The shift in the cricital fibre orientation 
angles is clearly shown in the Figs. 4.12- 4.20. Also it can be seen that the plots 
are no longer symmetric for plates with a cutout 
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Figure 4.6; Non-dimensionalized initial buckling load versus fibre orienta- 
tion (0), (Material- M2), a/b=:l. 



Figure 4.7: Non-dimensionalized initial buckling load versus fibre orienta- 
tion ( 0 ), (Material- M2), a/b=2. 
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Figure 4.8: Non-diiiieiisioiialized initial buckling load versus fibre orienta- 
tion (0), (Material- M2), a/b=3. 



Figure 4.9: Non-dimensionalized initial buckling load versus fibre orienta- 
tion ( 0 ) with a cutout , (Material- M2), a/b=l. 
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Figure 4.10: Non-dimensionalized initial buckling load versus fibre orien- 
tation ( 6 ) with a cutout , (Material- M2), a/b=:2. 



Figure 4.11: Non-dimensionalized initial buckling load versus fibre orien- 
tation { 6 ) with a cutout , (Material- M2), a/b=3. 
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Figure 4.12: Non-dimensionalized initial buckling load 
tation (B) with and without cutout , (Material- M2), 


versus fibre orien- 
a/b=l, a/h=10. 


70 


50 
40 
30 

20 

10 
0 

Figure 4.13: Non-dimensionalized initial buckling load versus fibre orieu 
tation ( 9 ) with and without cutout , (Material- M2), a/b=l, a/h=:50. 
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Figure 4.14: Non-dimensionalized initial buckling load versus fibre orien- 
tation {()) with and without cutout. , (Material- M2), a/b=l, a/h=100. 



Figure 4.15: Non-dimensionalized initial buckling load versus fibre orien- 
tation {()) with and without cutout , (Material- M2), a/b=2, a/h=:10. 
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Figure 4.16: Non-dimensionalized initial buckling load versus fibre orien- 
tation ((?) with and without cutout , (Material- M2), a/b=2, a/h=50. 



Figure 4.17: Non-dimensionalized initial buckling load versus fibre orien- 
tation ( 9 ) with and without cutout. , (Material- M2), a/b=2, a/h— 100. 
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Figure 4.18: Non-dimensionalized initial buckling load versus fibre orien- 
tation ( 9 ) with and without cutout , (Material- M2), a/b=3, a/h=10. 



Figure 4.19: Non-dimensionalized initial buckling load versus fibre orien- 
tation (^) with and without cutout , (Material- M2), a/b=3, a/h=50. 
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Figure 4.20: Non-dimensionalizcd initial buckling load versus fibre orien- 
tation {()) with and without cutout. , (Material- M2), a/b=3, a/h=100. 

4.9 Biaxial In-Plane Loadings and Ny 

This section deals with the study of stability envelops under biaxial loading on the 
response of laminated coinpositc plate with and without a cutout. The results for a 
plate without a cutout are compared with the results available in the literature. [13]. 
It is found that the results are in very good agreement with the results m the 
literature. 

Figs. 4.21-4.28 show the stability envelops under biaxial loadings for plate aspect 
ratio’s {a/b - 1,2) and Icngtli/thickness ratio’s 10 and 100. lU in the figures 
represent the noudimensionalised buckling load {lix — N^b'^/Eih ' ) and (Ry — 
NybyEth^). 

For a/b = 1 and a/h = 10 . Figs. 4.21 and 4.25 show the interaction curves 
without a cutout and with a cutout. The following observations can be made from 

the figures. 

• Except for cross ply laminates the stability envelops for all angle ply laminates 
are linear for a plate without a cutout. But this behaviour is not observed for 
the case of a plate with a cutout. 

. ^ = 45“ gives the highest value of (A/, and A/,) for any combination of the loads 
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in both the cases and for cross ply laminate {N^ and Ny) is the mininum. 

• Iherc is a sudden change in the slope ol stability envelops for almost all 
laminates for a plate with a cutout. Whereas when it comes to plate without 
a cutout, there is no such sudden change of slope . 

• From the stablity envelops shown in Fig. 4.21, it can be observed that for 8 = 
30° and 0 = G0° laminates there is a region between two points in which both 
the laminates have the same interaction curve. But when it comes to a plate 
with a cutout, 8 — 30° and 6 = 45° laminates has the same interaction curve 
at a single point. At this particular point, both these laminates withstand the 
same amount of buckling load. 

Figs. 4.22 and 4.26 show the interaction curves for a plate without a cutout and 
with a cutout. The following are the observations that can be made from the figures. 

• The stability envelops are linear or piece wise linear for a plate without a 
cutout whereas they are no longer linear for a plate with a cutout. 

• 0 = 45° gives the highest value of {Nj. and Ny) for any combination of the 
loads in both the cases and for cross ply laminate {NxandNy) is the mininum. 

• There is a sudden change in the slope of stability envelops for almost all 
laminates for a plate with a cutout. Whereas when it comes to plate without 
a cutout, the slope changes only for 8 = 30° and 8 = 60° laminate. 

• The portion for which [Nx + Ny) is same for both 8 = 30° and 8 = 60° 
laminates is less as compared to thin laminates. 

• Comparing the magnitude of the buckling loads for plate with and without 
cutout, we can observe that there is some 40 precentage reduction in the 
buckling load. 

For a/b = 2 and a/h = 10 (Fig.4.23 and Fig. 4.27) the following observations 
are made. 
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• Here for both the cases ie., plate without a cutout and with a cutout, the 
stability envelops are piecewise linear but for certain portions of the envelops 
there is a change in the slope. 

• The curves are almost perpendicular to the Y axis for one portion of them and 
are almost perpendicular to the X axis for the other, portion. This is due to 
the fact that at these portions of the curves, for a small decrease in the value 
of Nx the strength of the laminate to carry Ny increases rapidly. 

• From Fig. 4.22 we can observe that 0 = 45° and 6 = 60° laminates have 
interaction curve at one particular point. There is a reduction in the buckling 
load when it conies for a plate with a cutout. 

For a/b = 2 and a/h = 100 (Fig.4.24 and Fig. 4.28) the following observations 
are made. 

• Here for both the cases ic., plate without a cutout and with a cutout, the 
stability envelops are piecewise linear. 

• For different combinations of Nx and Ny , either 6 = 45° ov 6 = 60° laminate 
gives the maximum value of (A^x + Ny). 

• For crossply laminate, the slope of the stability envelop changes at high value 
of Ny for both the cases, ic., for plate without and with a cutout, for 9 = 45° 
laminate change in slope occurs at high value oi Nx whereas 6 = 60° laminate, 
change in slope occurs twice. One is when both Nx and Ny are high and 
another is when only Nx is high and Ny is small. 
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Figure 4.21: Compaxative study of stability envelop for biaxial in-plane 
compressive loads Nx and Ny, (Material- M2), a/b=l, a/h=10, without 
cutout, {b = 254mm, h = 2.1127nm). 



Figure 4.22: Comparative study of stability envelop for biaxial in-plane 
compressive loads Nx and Ny, (Material- M2), a/b=l, a/h=100, without 
cutout, (b — 254mm, h = 2.112mm). 
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Figure 4.23: Comparative study of stability envelop for biaxial in-plane 
compressive loads and Ny, (Material- M2), a/b=2, a/h=10, without 
cutout, {b = 254mm, h = 2.112mm). 



Figure 4.24; Comparative study of stability envelop for biaxial in-plane 
compressive loads Nx and Ny, (Material- M2), a/b=2, a/h=100, without 
cutout, {b = 254mm, h = 2.112m7n). 
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Figure 4.25: Comparative study of stability envelop for biaxial in-plane 
compressive loads and Ny, (Material- M2), a/b=l, a/h=10, with 
cutout, (b = 2547rim, h = 2.1l2rnrn). 



Figure 4.26: Comparative study of stability envelop for biaxial in-plane 
compressive loads Nx and Ny, (Material- M2), a/b=l, a/h=100, with 
cutout, {b = 254mm, h = 2.112mm,). 






Figure 4.27: Comparative study of stability envelop for biaxial in-plane 
compressive loads Nx and Ny, (Material- M2), a/b=2, a/h=10, with 
cutout, (6 = 254?nm, h = 2.112m7n). 



Figure 4.28: Comparative study of stability envelop for biaxial in-plane 
compressive loads Nx and Ny, (Material- M2), a/b— 2, a/h 100, with 
cutout, {b = 254mm, h = 2.1127mn). 





4.10 Biaxial In-Plane Loadings and N^y 

III this section stabily envelops a.re shown for both crossply and aiigleply symmetric 
laminates [±/?] 2 .s with plate aspect ratio’s 1 and 2 and length/thickness ratio’s 100 
s-iid 10 for both scpiaie and rectangular laminates. Figures 4,29 "4.39 shows the 
interaction curves for and N^y 

For a/b=l and a/h^-^lOO (Figs. 4.29 and 4.31) the following observations are 
made. 

• Stability envelops are not symmetric about N^y = 0 line except for a crossply 
laminate for both the cases, ie.,with and without a cutout. It is observed that 
the value of {N^ + \ Nxy\) for initial buckling load is more for positive shear as 
comapared to negative shear. This trend is also observed in the plate with a 
cutout. 

• Stability envelops 6 = 30° and 9 = 60° laminates are comparable though 
9 — 30° laminate gives slightly higher values of {Nx + jA^iyj)- 

• From the figures we could observe that there is no such unique point at which 
two or more laminates carry the same amout of axial compressive and shear 
loads. 

For a/b=l and a/h=100 (Figs. 4.30 and 4.32) the following observations are 
made. 


• Stability envelops are not symmetric about Nxy — 0 line except for a crossply 
laminate for both the cases, ie.,with and without a cutout. It is observed 
that the value of {Nx + |iV,,„|) for initial buckling is more for positive shear as 
comapared to negative shear. This trend is also observed in the plate with a 
cutout. 

• For negative shear loading the stability envelops of 9 = 30° and 9 = 45° 
laminates are comparable with each other while 9 = 60° and crossply laminates 
are comparable with each other. The same trend can be observed in a plate 
with a cutout. 
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• For both the cases the interaction curves for the laminates intersect each others 
at some unicjne points. At these points dillerent laminates carry the same 
amount of biaxial loads, and N„j. In other words it can be said that there 
are some unicpie combinations of in-plane compressive and shear loadings that 
can be carried out by different laminates. 

• From the figures we cmild observe that there is a reduction in the buckling 
load in all the laminates when a cutout is placed in the plate. The variation 
varies for different laminates. 



Figure 4.29: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads iV® and Nxy, (Material- M2), a/b=l, a/h=10, 
without cutout, (b = 254rnm, h = 2.1127nm). 
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Figure 4.30; Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads and N^y, (Material- M2), a/b=l, 
a/h=100, without cutout, {b = 254mm, /i = 2.112mm). 



Figure 4.31: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads Nx and Nxy, (Material- M2), a/b=l, a/h=10, 
with cutout, {b = 254m7n, h = 2.112mm). 
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Figure 4.34: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads Nx and Nxy with and without cutout, 
(Material- M2), a/b=l a/h=10, (6 = 254mm, = 2.112mm). 



Figure 4.35: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads Nx and Nxy with and without cutout, 
(Material- M2), a/b=l a/h=10, {b = 254mm, /i = 2.112mm). 
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Figure 4.36: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads and N^y with and without cutout, 
(Material- M2), a/b=l a/h=100, (6 = 2547nm, h = 2.112mm). 



Figure 4.37: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads and N^y with and without cutout, 
(Material- M2), a/b=l a/h=:100, {b = 254mm, h = 2.112mm). 
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Figure 4.38: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads Nx and Nxy with and without cutout, 
(Material- M2), a/b=l a/h=100, (6= 254mm, /i = 2.112mm). 



Figure 4.39: Comparative study of stability envelop for biaxial in-plane 
compressive and shear loads Nx and Nxy with and without cutout, 
(Material- M2), a/b=l a/h=100, {b = 25477im, h = 2.112mm). 
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4.11 Effect of Cutout to Plate Ratio on the Buckling 
Load of a Laminated Composite Plate. 

In this section the efiect of cutout to plate ratio on the buckling load of two laminates 
([0/90)5 and [—45/45)5) is made. The material properties used are M-2 and the 
results are plotted for in-plane compressive loading N^. The plot for [0/90)5 laminate 
is shown in Fig. 4.40 and that for [-45/45)5 is shown in Fig. 4.41. The figures show 
the variation of the buckling load for various hole to plate ratio’s in X direction 
with the hole to plate ratio in Y direction. Also the corresponding buckling load for 
a plate without a cutout is shown in the plots. The following observations can be 
made from the figures . 

• For [0/90)5 laminate the buckling load increases gradually until the hole to 
plate ratio in Y direction reaches 0.3 for c/a — 0.1 and there after it drops 
down. But for a [45/ — 45)5 laminate the maximum value of the buckling load 
occurs when the hole to plate ratio in Y direction reaches 0.38. 

• Considering the case when c/a = 0.2, its interesting to note that the buckling 
load exceeds the value without a cutout. The same trend can also be observed 
with [45/ - 45)5 laminate. Here the buckling load reaches it peak with d/b = 
0.18 and there after it comes down. So, for this hole to plate ratio the stiffness 
of the plate is increasing than that of a plate without a cutout. 

• For the remaining hole to plate ratio’s, the buckling load gradually drops down 
for both the laminates. 




Figure 4.40; Buckling load vs Hole to plate ratio in Y-direction 
in-plane load , (Material- M2), a/b=l a/h=:10 



Figure 4.41; Buckling load vs Hole to plate ratio in Y-direction 
in-plane load Nx , (Material- M2), a/b=l a/h=10 


under 


under 
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Chapter 5 


Conclusions and Future Scope 

5.1 Conclusions 

In the present study the initial buckling load of the composite laminates with and 
without a cutout is studied for different fibre orientation angles, plate aspect ra- 
tios and length/thickness ratios. The computations are done for both uniaxial and 
biaxial loads (both compressive and shear). Interaction curves under the above 
said loadings are made taking the efl'ects of the cutout. Finite element method is 
employed to compute the eigen values which are obtained from the finite element 
solution. The fundamental eigen value corresponds to the initial buckling load of 
the laminated composite plate. 

The following conclusions are made from the present investigation. 

• The effect of transverse shear is negligible for thin laminates, but coming to 
thick plates it has to be duly taken in to account for calculating the buckling 
load. Hence HSDT should be used for moderately thick and thick plates. 

• The effect of laminate orientation angle on the initial buckling load of an 
antisymmetric laminate decreases with the increase in the length/thickness 
ratio. 

• For thick laminates the eflFect of laminar orientation angle of antisymmetric 
laminate on the initial buckling load decreases as the aspect ratio increases. 
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• For thick larninatos the critical fibre orientation angle becomes smaller as the 
aspect ratio of the laminate increases. 

• The buckling load of a laminated plate decreases when a cutout is placed in it 
at the center. There is a shift in the critical fibre orientation angle for a plate 
with a cutout. 

• Interaction curves for and for a plate without a cutout are linear or 
piecewise linear. But with a cutout, the linear relationship no longer exists. 

• There are some unique points where the laminates carry same amount of buck- 
ling load. This combination exists even for a plate with a cutout. But the 
magnitude of buckling load is considerably reduced when a cutout is placed. 

• The stability envelops for combined shear and compressive loadings ( and 
Nxy ) are symmetric only for a crossly laminate. But for other laminates they 
are no longer symmetric. The same trend can also be observed for a plate with 
a cutout . 

• There is a decrease in the buckling load of a composite plate when a cutout 
is placed. But there exists one specific hole to plate ratio where the buckling 
load exceeds that of a plate without a cutout. So we can say that the stiffness 
of the plate increases with the material removal. 

5.2 Scope for Future Work 

The areas which require further study in the future are given below. 

• The analysis of plates with multiple cutouts of different geometries and at 
different locations can be made. 

• The laminate orientation angles can be optimized to get the optimistic buckling 
load. 

• Response of laminates with plydrops need to be investigated. 
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• The present analysis is done taking only linear relationships. So further inves- 
tigation can be made taking in to effect the non-linear effects. 

• Stability analysis of thick laminates with different combinations of thermal/hygroscopic 
and static in-plane loading. 
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APPENDIX A 
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in which, 

= = 1 , 2 , 3 , 4 , 5, 7 

and i/s = [{zfc — Zk-i} — ^{z^k — -z^-i} + — ^°fc-i}]. 
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The aforcmcnlioiicd shape fimclious can be expressed in natural co-ordinates as; 

- 1 )M-? - 2 )(i) - 1 )^-, - 2 ) 

^2 = ^ (f + U’e? - 2 )(>; - - 2 ) 

W' = i({ -b !)=(? - 2),(,, -b 1)*(, _ 2) 

W; = i(? - - 2)(,, + 1)^{, - 2) 

w; = i({ - l)*(-n - 2 ) 

W; = - -!)(>?- 1)^-1 - 2 ) 

w; = - ^(« + l)ne -!)(>) + - 2 ) 

JV; = - l?{-( -!)(>)+ Ifil - 2) 

Wj = i(« - l)"(-f - 2){, - 1)2(-.7 - 1) 

N'la = ^(« + !)"(« - 2)(’l - !)’(->) - 1) 

W;, = - 2)(.) -b l)"(.l - 1) 

Wjj = - - lfl-( - 2) (r, -b Ifiv - 1) 

/v ;3 = - 1 )^-? - ixi - - 1 ) 

w;, = -b !)’(« - i)(’i - ijM-i - 1) 

= ^(« + !)"(« - M*’' + 

ivie = - -!)"(-«- 1)(>J +!)"(’)- 1) 
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APPENDIX C 

The details of [Bq] and [Bg] matrices are as shown: 
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In the above matrices, for columns 1 and 2, i = 1 — 12 (shape functions cor- 
responding to uq and vq) and for columns 3 and 4, i = 1 — 16 (shape functions 
corresponding to Wb and Ws)- 
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